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We study dynamical effects in proton breakup from a weakly bound state in an exotic nucleus 
on a heavy target. The Coulomb interactions between the proton and the core and the proton 
and the target are treated to all orders, including also the full multipole expansion of the Coulomb 
potential. The dynamics of proton nuclear and Coulomb breakup is compared to that of an equiva- 
lent neutron of larger binding energy in order to elucidate the differences with the well understood 
neutron breakup mechanism. A number of experimentally measurable observables such as parallel 
momentum distributions, proton angular distributions and total breakup cross sections are calcu- 
lated. With respect to nuclear breakup it is found that a proton behaves exactly as a neutron of 
larger binding energy. The extra "effective energy" is due to the combined core-target Coulomb 
barrier. In Coulomb breakup we distinguish the effect of the core-target Coulomb potential (called 
recoil effect), with respect to which the proton behaves again as a more bound neutron, from the 
direct proton-target Coulomb potential. The latter gives cross sections about an order of magnitude 
larger than the recoil term. The two effects give rise to complicated interferences in the parallel mo- 
mentum distributions. They are instead easily separable in the proton angular distributions which 
are therefore suggested as a very useful observable for future experimental studies. 
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I. INTRODUCTION 



'T^ In a recent paper Liang et al.llj studied experimentally dynamical effects in the Coulomb breakup of ^^F. The 

^ motivation was that Coulomb dissociation, being the inverse of the radiative capture reaction, is a useful technique 
^ for studying stellar nucleosynthesis involving short-lived nuclei where direct measurements are difficult [7 . The authors 

'~~' argued that with an increasing number of radioactive isotope beams available for studying capture reactions that occur 
, in stellar environments, such as those in the rp-process, more measurements will be performed by Coulomb dissociation 

^ because of the short lifetimes that make direct measurements impractical. In order to extract reliable information 
^^ on radiative capture the reaction mechanism in Coulomb dissociation has to be understood in detail. Furthermore 

\^ Liang et al. [1] stressed the importance of the dynamic polarization effect which they interpret as a displacement of 
CN the valence proton behind the nuclear core and a subsequent shielding from the target. This effect manifests itself 

^~l as a reduction of breakup probability compared to first-order perturbation theory predictions. Coulomb dissociation 
f — . experiments [3]-[5] have been studied with different approaches, among which first-order perturbation theory is indeed 
^^ often used with final plane wave functions. While Coulomb dissociation of loosely bound neutron-rich nuclei has been 
'""' studied extensively and is fairly well understood [51 [7] , this is not true for proton-rich nuclei in which the loosely 

\ I bound valence protons actively participate in the reaction. It has also been suggested that the inclusion of higher- 
order corrections is required [HI [9] . Essentially exact calculations of realistic three-body models of break-up based on 
the Faddeev equations including the Coulomb potential exist for few body nuclei [TU] . These methods have also been 
applied to neutron break up from medium mass exotic nuclei using cluster models |llj and it is expected that they 
^ should perform equally well for proton breakup. 

Besides the well known astrophysical implications proton rich nuclei present a number of unusual features such 
as the two proton radioactivity or beta-delayed proton emission which make them very appealing to study and to 
compare to the neutron rich nuclei. An account of the richness of the "Physics of the proton rich side of the nuclear 
chart" can be found in Ref. 12 . 

In Ref.[T] ^^F + ^^Ni and ^^F + ^'^'^Pb breakup angular distributions measured at lOA.MeV were discussed and 
compared to first-order perturbation theory prediction for Coulomb breakup and to fully dynamical calculations. At 
the end of their paper the authors used first-order perturbation theory to calculate proton breakup via an increase 
of the effective binding energy of the valence proton treated as if it were a neutron, according to a model that was 
proposed some time ago by some of us [13]. Liang et al. T; found that our idea did indeed work and the angular 
distributions could successfully be reproduced using an increase in the binding energy of 1.2 MeV, contrary to our 
model suggestion of 3.2 MeV. 

We were intrigued and interested by such a result: first we were happy to notice that our very schematic model 
would work in comparison with real data; second we decided that it was worth studying if and how the model could 
be made more quantitatively reliable. This paper reports results of such a study for some observables such as parallel 



momentum distribution and proton angular distributions following breakup. A forthcoming paper will discuss core 
angular distributions. 

In the next section we recall the model for nuclear and Coulomb breakup for a valence nucleon treated to all orders 
introduced in 14J, that we developed successively in Ref.[15_ to treat protons. Our starting point is, as in all our 
previous works, a first order perturbation theory amplitude such as that given by Eq.(l) of J16j . In that paper the 
amplitude was used for transfer to bound states, thus the neutron-target potential acted only to first order. On the 
other hand when the same formula is used for transitions to final unbound states |17j . in order to describe breakup, 
the final state is a proper scattering state with respect to the target and thus all re-scattering terms are included 
[inil^n]- We then compare in Sec. Ill results of the full calculation to those of perturbation theory. The amount and 
characteristics of the breakup due to the core-target and to the proton-target Coulomb repulsions are then studied 
separately. We do that at fixed impact parameter values in order to understand the reduction of breakup probability 
for protons as compared to that of neutrons of the same binding. This is because our model predicted that the 
effective binding energy would depend on the distance between the two nuclei, being due to the mutual projectile- 
target Coulomb barrier. The impact parameter dependence enters also in the model of Ref. |15j where we showed 
that the Coulomb breakup could be treated by time dependent perturbation theory at large impact parameters, while 
it could be described by an all order eikonal model at small impact parameters. Sections III contains results for 
®B, while Section IV is devoted to ^^F. In Section V the proton and neutron angular distributions due to Coulomb 
breakup are shortly described. Finally our conclusions are given in Sec. VI. 



II. FORMALISM 




FIG. 1: (Color online). Coordinate system. 



The Coulomb potential responsible for proton breakup is 
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where Vc = Z^Zte^, Vy = ZyZfC^ and Vq = {Zy + Z^Z^e^ . /3i and /32 are the mass ratios of proton and core, 
respectively, to that of the projectile. Zc, Zj and Z„ are the core, target and proton charge respectively. The 
coordinate system used in this paper is shown in Fig. IT] Following [TH [TJ] , the perturbation theory phase is defined 
as Xvert = \ /<iie*"*y(r, t), which calculated explicitly gives 
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with by = bc + r±, to = {ef — Eq) /h and Eq is the neutron initial bound state energy while £/ is the final neutron-core 
continuum energy. Since Vq — Vc + Vv, Eq. ([2]) can be written as 

x'^xWi,V,) + x{-l32,V,) (3) 

where each term is given by 

XW, V)^^ [e^^'^^/^Koitub/v) - KoiuRjv)) , (4) 

x(/3i, Vc) describes the recoil of the core whereas x(~/32, Vy) accounts for the direct proton-target Coulomb interaction. 
The latter vanishes in the case of the neutron. 

The expansion of Eq. (HI) to first order in f yields the dipole approximation to the phase: 

X^ - z {Kq{uRj_/v) 
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The constant factor is now (K/Si — K/32) while for a neutron would be Vb/3i [T5] . 

We recall here that in Ref.15] it was suggested that proton breakup could be treated as neutron breakup and 
the Coulomb potential of Eq. (fl]) could be approximated by a dipole expansion if the neutron was given an effective 
binding energy which would take into account the combined effect of the core and target Coulomb barriers. As we 
mention in the Introduction, this model has been applied in Ref. [T] . Liang et al. used first-order perturbation theory 
to calculate proton breakup via an increase of the effective binding energy of the valence proton treated as if it were 
a neutron. The authors fitted their the angular distributions using an increase in the binding energy of 1.2 MeV. We 
shall show in the following that accurate calculations of proton vs neutron breakup suggest indeed that in certain 
cases proton breakup can be understood in terms of a more bound neutron breakup. However we believe now that a 
better estimate of the effective energy is given by 

which through the factors /3i and /32 is consistent with Eq.(fl]). R^ the position of the projectile top of the Coulomb 
barrier, d is the distance between the center of the two nuclei for which the top of the two Coulomb barriers of 
projectile and target coincide. 

Our expression for the differential cross-section is 

where \Sct{bc)\'^ is the core-target elastic scattering probability. 

The probability amplitude in Eq.([7]) has been written as the sum of three pieces: the recoil term, 

g-= = J dre'^'^-^U^i^^ i°s ^ _ 1 _ ^fl£ log |^ 

+ix(/3i,K)), (8) 

obtained in the sudden limit according to P^ [T5] in order to include all orders in the final state Coulomb interaction 
of the core with the target. Similarly, the second term in our probability amplitude is the direct proton Coulomb 
interaction. It has the same form as Eq.Q but for the substitution Vc — >■ K, ^c — >■ by and /3i — )• — /32. Both the 
direct and recoil term contain a regularization of the first order term divergency consisting in substituting it with the 
corresponding first order perturbation theory term. Such regularization method was first proposed in Ref. |20j and 
used also in Refs.^HHl]. 
Finally, the nuclear part is 

guuc ^ I rffe-*^>,(7=0 (e^X"'*'''') - l) . (9) 



in terms of the one-dimensional Fourier transform in z— direction 
gets: 



If the Coulomb part of the amplitude g'~^°'^ = g'^'^'^ + g''-'^^ is simply expanded to first order in x and it is written, 

then making also the dipole approximation one 
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leading to symmetric momentum distributions. 

A number of experimental data for neutron breakup have been analyzed by separating high impact parameters 
such that perturbation theory would be valid. We shall try in the following to do the same for proton breakup and 
in order to clarify the effects of the Coulomb interaction we will show separately results for the Coulomb or nuclear 
probability only. Furthermore we will show momentum probability distributions at fixed impact parameters, namely: 



dPjbc) ^ 1 
dk^ 87r3 
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|(7(&c)P will be the nuclear or Coulomb amplitude and in the latter case we will indicate in which approximation the 
calculation was done. 



III. APPLICATION TO *B 



8B+208Pb at 72A MeV in Coul. Breakup reaction to aii order 
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(Color online). Recoil term, Eq. 



red dashed line and direct term, green dotted line. The black full line contains 



In our previous work in which the method presented here was introduced, we applied it to study a reaction at 
936A.MeV, the energy was high enough to justify the eikonal approximation. Before starting the calculations presented 
in this and the following section we have compared our results for the total breakup probabilities against the values 
provided by Ref . [18] . We have found that for nuclear breakup our results agree with the dynamical calculations down 
to lOA.MeV if we substitute the impact parameter by the distance of closest approach of the corresponding classical 
trajectory. For Coulomb breakup our probabilities decrease with the impact parameter as those of Ref.|18j but they 
are about 50% larger. This is due to the use of a straight line relative motion trajectory in our method. For nuclear 



TABLE I: Barrier radii, initial binding energies and effective energy parameters for a ^"^^Pb target. 
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breakup the substitution of the impact parameter with the distance of closest approach is enough to cure the method, 
because the nuclear form factors are different from zero only in a small region around the distance of closest approach 
where the classical trajectory is well approximated by a straight line. On the other hand the same cannot be said 
of the Coulomb form factor. We have checked however that our method starts to give reliable absolute probability 
values around 70-lOOA.MeV which agrees with what we found in our previous works on neutron breakup P^I20) . 
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FIG. 3: (Color online). Proton vs neutron wave functions for a P3/2 single particle state in a *B and d5/2 in ^^F as indicated. 



Thus the formalism described in the previous section has been applied to proton breakup of ^B and ^^F against 
Pb target at a beam energy of 72A.MeV which is a typical energy used in several laboratories and for which our 
results should be reliable. The projectiles are taken as a two-body object. Radial wave functions have been obtained 
by numerical solution of the Schrodinger equation in Woods-Saxon potentials with depths adjusted to reproduce the 
experimental separation energies (0.137MeV and 0.6MeV respectively). The radius parameter of these Woods-Saxon 
potentials has been taken as 1.3 fm and the diffuseness as 0.6 fm. 

In Fig(2] we show separately the effect of the recoil, Eq. (l8| term, red dashed line while the green dotted line 
represents the direct term obtained from g'^^^ . The black full line contains both. It is clear that the interference 
between the two is mostly destructive and it becomes constructive only at very large impact parameters. This effect 
is the equivalent in our model of what Liang et al. [1] call the "shielding effect" of the proton by the core. The direct 
term alone, being proportional to (^2 is indeed in absolute value always larger than the recoil term. However it is the 
effect of the interference between direct and recoil terms that causes the reduction of the Coulonrb breakup in the 
proton halo case. 

Following Refs.p and [13] we have then calculated the Coulomb breakup to all orders (regularized) and in the 
dipole approximation for a neutron with several different binding energies. Our hypothesis, that one could use a 
neutron wave function corresponding to an effective separation energy Sn larger than that of the proton, has been 
checked for various values of 5„ and we have found that for Sn = 0.6 MeV the model works quite well. This is a factor 
four smaller than the value Sn — 2.4 MeV predicted in [13]. Liang et al. in [1] found indeed that in the description 
of the breakup of ^^F the dipole approximation and perturbation theory to first order could be used by introducing 
a neutron wave function corresponding to an effective energy smaller by about a factor three than the one of our 



8B+208Pb at 72A MeV proton vs neutron parellel momentum distribution 
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FIG. 4: (Color online). Differential cross sections for neutrons of different separation energies and proton breakup cross section, 
all normalized to one. The integrated values are given in the legend. 



original model. The reason is that the model of [IB] used an effective binding energy given by an intuitive but not 
very accurate expression. In this work we use instead Eq.(l6]) which through the factors /3i and /32 is more consistent 
with the Coulomb potential definition. In order to understand the proton vs neutron breakup dynamics we start by 
showing in FiglslLIIS the single particle wave function {Sp = 0.137 MeV) for a P3/2 proton with respect to a ^Be core 
by the full blacT line, for a neutron with the same Sn by the red dashed line and finally for a neutron with Sn — 0.6 
MeV by the green dotted line. The RHS is for ^^F wave function. 

We have made all order proton calculations and calculations for a neutron of separation energy Sn = 0.137 MeV and 
Sn = 0.6 MeV. In the latter case we have studied the impact parameter dependence of the momentum distributions 
obtained to all order and compared the results to those corresponding to the dipole approximation. It is interesting to 
note that as discussed in |13) , and in particular in the appendix of that work, the dipole approximation for a neutron 
of large binding energy reproduces quite well the proton exact calculation and becomes actually identical to it at large 
impact parameters. 

At this stage we remind the reader that in the case of breakup of a neutron both the " small" width of the momentum 
distribution and large breakup cross sections were due to the small separation energy and corresponding long tail of 
the wave function. The proton breakup characteristics are instead quite different as we can interfere from FigH] It 
shows the proton breakup parallel momentum distribution from ^B on a ^''^Pb target compared to three distributions 
from an hypothetical neutron bound by three possible energies as indicated. The distributions are normalized to one. 
The total breakup cross sections including nuclear and Coulomb interactions to all orders are given in the legend. 
The smallest neutron separation energy case provides the closest cross section to the proton's one while it is in the 
largest separation energy case that the neutron distribution width gets closest to the proton distribution width. This 
is consistent with the fact that it is the tail of neutron wave function corresponding to S„=0.6MeV in Fig|3] which is 
closest to the proton's wave function tail. 

Therefore we understand that as in the neutron case, the large cross sections are associated with small separation 
energy but the corresponding momentum distributions have large width because the wave functions do not have an 
extended tail. The model of a neutron of larger separation energy representing a proton of smaller separation energy 
works therefore well for nuclear breakup. In the case of the proton Coulomb breakup the results are not easy to 
interpret instead because direct and recoil terms interfer mostly destructively (cf. Fig(2]). However the interference 
becomes constructive at the very large impact parameters where the recoil term is negligible. There, a large breakup 
can be seen as a consequence of the direct term which as Eq.(l5]) shows, being proportional in first order to ^2, will be 
large if the proton is well displaced from the core. 

We have also studied the spectra corresponding to the m=0,l components of the wave function for the direct 
and recoil terms in Coulomb breakup including interference and we compared to the same quantities in the case 
of a neutron. The m=0 componenent gives asymmetric distributions while the m=l component gives symmetric 
distributions. These results indicate that small polarization effects can lead to asymmetric spectra in the proton case 
and that this is due to the large direct breakup term. Neutron spectra do not show this effects. 



IV. APPLICATION TO ^^F 



We start this section by suggesting the reader to look at the proton wave function for the d5/2 ground state of ^^F 
back in Fig. [3] . The same figure contains two neutron wave functions for separation energies as indicated that we 
have used in the calculations shown in the subsequent FigsjSJIHlandlT]. 



1 7F+208Pb at 72 A MeV stripping reaction 1 7F+208Pb at 72 A MeV diffraction reaction 
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FIG. 5: (Color online). Diffraction and stripping terms of the nuclear breakup of the proton d5/2 ground state in ^^F 

By looking at the wave functions in Fig. |3] and the nuclear breakup diffraction and stripping terms, Figjs] we see 
that the "neutron-like" model works well also for the d5/2 ground state in ^^F at 72A.MeV incident energy. The best 
"model" separation energy here seems to be 1.7MeV. This value is larger than what suggested in [T] but smaller than 
the value predicted in P^. It agrees well with our new effective energy estimate Eq.Q. 



17F+208Pb at 72 A IVIeV proton vs neutron parellel momentum distnbution 
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FIG. 6: (Color online). Differential breakup cross sections for neutrons of different separation energies and the proton (d5/2 
ground state in ^^F) cross section. All normalized to one. The integrated values are given in the legend. 
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FIG. 7: (Color online). Direct and recoil terms in Coulomb breakup for the m=l component of the ^^F wave function. 

To finish tiiis section we show in FiglHlthe total cross section distribution for the breakup of the A^/2 proton ground 
state in -'^^F at 72 A. MeV, compared to two neutron breakup distributions. What is evident here is a strong asymmetry 
in the spectrum which instead does not appear in the neutron case. Also in this case the presence of a direct Coulomb 
breakup term is extremely important because the recoil term is strongly reduced by the centrifugal barrier of the d- 
state. The immediate consequence is that now the interference between direct and recoil term is mostly constructive. 
The asymmetry in the spectrum originates from the interference of direct and recoil term which is present in the low 
(m=0,l) components of the initial wave function which dominate the breakup spectrum. Such effects are shown by 
Fig. [7]which contains the probability spectra for the m=l projection of the initial angular momentum wave function. 
The m=2 component gives a symmetrical distribution not shown here. 



ANGULAR DISTRIBUTIONS 
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FIG. 8: (Color online). Proton and neutron angular distributions after breakup. Details are given in the legend. 



9 

This section is devoted to the discussion of the angular distributions of the breakup protons presented in Fig|8] 
Here again we show the total result for protons and at the same time the direct and recoil terms. We show also 
two calculations for neutrons. One for neutron having the same binding energy as the proton and the other for the 
effective energy previously obtained. It is very interesting to see that the recoil term of the proton calculation and the 
calculation for a neutron with the effective binding energy give very close results. This confirms the interpretation that 
the core-target Coulomb effect is basically a recoil effect for which the difference between a neutron and a proton is that 
the proton looks "more effectively" bound. On the other hand the direct term gives the most important contribution 
to the proton breakup, its interference with the recoil term is sometimes constructive, sometimes destructive. This 
term is not present in the neutron case and it is at the origin of the special dynamical effects of proton breakup 

Finally we notice that a neutron of the same binding energy as the proton would give quite a large breakup in the 
case of the p-state of ®B, thus suggesting a halo-like behavior for this nucleus, while for ^^F a neutron of the same 
binding would anyhow give much less breakup at most angles. We can understand this effect by looking at the wave 
functions in Figlsl The tail of a neutron wave function for ®B, being a p-state, is much more pronounced than the tail 
of a possible d-neutron wave function in ^^F. Therefore there is no "structural" halo effect in the Fluorine case but 
the large proton breakup probability in this case has a dynamical origin being due to the strong direct proton-target 
repulsion. 

VI. CONCLUSIONS 

In this paper we have calculated proton and neutron Coulomb and nuclear breakup from ^B and ^^F at 72 A.MeV on 
a 208p|-, target. Our method is an all order formalism based on the eikonal approximation, with a regularized first order 
Coulomb term. We have compared proton to neutron parallel momentum distributions and angular distributions. 
The neutron cases have been calculated for an "hypothetical" neutron having the same separation energy and angular 
momentum as the proton in the two projectiles and compared to calculations for a neutron having an effective binding 
energy larger than the true proton experimental binding energy. The effective values take into account the combined 
effects of the projectile and target Coulomb barrier. We have given an explicit formula to evaluate the effective 
energy which is in agreement with phenomenological findings. The model has been used to understand the origin 
of the strong dynamical effects seen in Coulomb breakup data. Our results clearly show that as far as the nuclear 
breakup mechanism is concerned the proton behaves as a neutron of larger separation energy. On the other hand 
parallel momentum distributions are affected by the recoil due to the neutron-core Coulomb potential and by the 
direct repulsion due to proton-target Coulomb potential, in a complicated way. This happens both in the calculations, 
as well as in the data. There are interference effects between the direct and recoil term which depend on the impact 
parameter. Thus even if they can be understood in the theoretical calculation, they would be difficult to disentangle 
in an experiment. Effects on the core angular distributions are also complicated and will be discussed in a forthcoming 
publication. However we have finally shown that the proton angular distributions would be a very good observable 
to study experimentally in order to separate the two Coulomb effects. Our results indicate that the proton angular 
distributions are dominated by the direct term while the recoil term gives almost an order of magnitude smaller cross 
section as a function of the proton angle. A neutron of the appropriate effective binding energy would have exactly 
the same angular distribution as that given by the recoil term for the proton. On the proton angular distributions 
the effects of interference seem to be negligible and we have found some only in the very forward angle region for the 
^^F projectile. We have also concluded that a large breakup cross section for a weakly bound proton can originate 
from dynamical effects without corresponding to an unusual structure of the bound state. Thus we have clarified the 
way in which the Coulomb barrier affects weakly bound protons in exotic proton rich nuclei off the stability line. We 
hope that our results would give some guideline not only in planning future breakup experiments for spectroscopic 
and astrophysical studies but more in general in studying the interplay between structure and dynamics of nuclei with 
weakly bound protons. 
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